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Motivation: Monte Carlo simulation and Sign problem 


Yuya Tanizaki (RIKEN BNL) Fermionic sign problem 


Exponential complexity to solve QFTs 


Quantum systems: the Hilbert space H and the Hamiltonian Å. 
Accept the lattice model of QFTs, then 


dimH = a. Volume i 


= Solving QFT completely (= diagonalization of H) is exponentially 
hard. 


Notice: We do not always need to know about all the states precisely. 
We are usually interested in 
Vacuum and several excitations, Thermal states, etc. 
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Path integral and Monte Carlo simulation 
To characterize the thermal state, we can use path integral, 
Z - try End = | DAexp(-5[4)). 
M 
The integration domain M (of an SU(N) gauge theory) is 
M= SU(N) Volume, 


We want a tool to evaluate Z without exponential complexity. 


Monte Carlo method: Consider the case e7314] > 0. 
Generate an ensemble {A;}; following e314, and evaluate 
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More about Monte Carlo simulation 


Do we really circumvent exponential complexity using MC method? 


- Typical values of |O(A;)| 
a l 
It indeed solves the exponential complexity for operators satisfying 


(O(A)) 
Typical values of |O(A;)| 


N 
1 
— > A; 
Error of N2 O(A:) 


~ (B - Volume). 


It has been quite successful to understand Hadron structures, 
thermodynamics of finite-temperature QCD, etc. 
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Sign problem 


What happens if STA) takes complex values so that e7514] 4 0? 


Reweighting: To use Monte Carlo method, we generate the ensemble 
LAj Ji following the phase-quenched distribution e~Re(SI4); 


(OA) oa. a. Pia Ase 
Cee) a > + = ae e-ilm(5[4]) 


(O(A)) = 


Since SIA) œ 8 : Volume, (e~!(SIAD), | = e76 Volume Af. Thus, 
Necessary N of configurations > e2? Volume As, 


Exponential complexity revives due to the sign problem. 
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Example of sign problem: Finite-density QCD 


Partition function of QCD is 


Zaco(T.m) = | DA Det(P(A, pa) + m) exp (Sua). 


quark gluon 


Sign problem: Det(D(A, py) +m) #0 at uq £ 0. 
Neutron star 
ə Cold and dense nuclear matter 
@ 2Msun neutron star (2010) 


ə Gravitational-wave observations 


Binary NS merger (NASA) 
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Method: Path integral on Lefschetz thimbles 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) - [Zoi (= + ar) 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a= 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (0 = 1,2) for the Airy 


integral: 
Ai(a) Yon f “an a 
= ya — —+az). 
= J, 2T p 3 


Ng: intersection number of the steepest ascent contour Ke and R. 
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Multiple integrals on Lefschetz thimbles 


Oscillatory integrals with many variables can be evaluated using the 


“steepest descent” cycles Js: (classical eom S'(z0) = 0) 


f d"™x eS) — wam | dee"), 
M 


o o 


Unlike one-dimensional case, the steepest descent manifold is not 
uniquely defined. 
=> Use of the homology H„ (Mc, {eP < 1}) becomes quite 


essential: 
> ZIJ), 


XOZ [Ka]. 


K 


H,,(Mc, (ep Ret) < 1} 


Hp (Me (e0 YO < 1}) 


K 


[Pham, 1967; Kaminski, 1994; Howls, 1997] 
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Multiple integrals on Lefschetz thimbles 


Concrete construction Pick up the Kahler metric ds? = ggdz’ @ dzi, 
and consider the gradient flow: 


Jo are called Lefschetz thimbles, and Im[S] is constant on it: 


In = Jo) 


im 2(t) = to) 


t--00 


Similarly, Ko = {z(0)|z(oo) = ze}. 
[Pham, 1967; Kaminski, 1994; Howls, 1997, Witten, arXiv:1001.2933, 1009.6032] 
[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
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Monte Carlo algorithm on Lefschetz thimbles 


Monte Carlo algorithms on J, has been developing since 2012. 


e Langevin on J, Cristoforetti, Di Renzo, Scorzato, 1205.3996 

ə Hybrid MC on J> Fujii, Honda, Kato, Kikukawa, Komatsu, 
Sano, 1309.4371 

ə Contraction algorithm Alexandru, Basar, Bedaque, 1510.03258 


These methods generate ensembles {z;}; on J, following e~ 5: 


d”zi 
' x, BOC) 
= d"zO(z) exp(-h715(2)) ~ N DDM 1 Ja” z,] i 
Z YA N LY 1 dna] 


Comment: la satisfies resurgence (Berry, Howls, '91, Howls '97). 
This numerical technique will allow to access analytic properties in 
the Borel plane for QFTs. 
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Lefschetz thimbles with fermionic determinant 
We would like to apply Lefschetz thimbles to treat the sign problem of 
Z= f DA e- Svmtar Ha Det (P(A.tig) +m), 

M 


Because of fermionic determinant, J, have their boundaries at 
Det(D(A, ug) +m) =0, or Re(Sym[A]) = too. 
The formula, [M] = $`, (M, K,)[J], still holds (1412.1891, 1412.2802 with 


Kanazawa). 
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Applications 


Applications: Case studies of fermionic sign problem 
e@ One-site Hubbard model (1509.07146, with Tomoya Hayata, Yoshimasa Hidaka) 


ə Multi-flavor massless QED» (to appear, with Motoi Tachibana) 
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AIET One-site Hubbard model 


Path integral for one-site Hubbard model 
We consider the (0 + 1)-dimensional fermion model, 


s= far (0/2 u= iole). 


The path-integral expression for the one-site Hubbard model 


(2 = 3 fy ard(r)): 


p f i 2 B2 
Zan | do (1 + chip t+ ¥/2))? 9897/20, 
2m7U R AB 


Fermion Det 


Integrand has complex phases causing the sign problem. 


p is an auxiliary field for the fermion number density: 
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AJETET One-site Hubbard model 


Behaviors of number density and Lefschetz thimbles 
Number density n with exact result and one-thimble result: 


-0.5 0.5 1 15 
(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


> 
JU 
Fujii, Kamata, Kikukawa,1509.08176, 1509.09141, Alexandru, Basar, Bedaque,1510.03258.) 


Lefschetz thimbles with —0.5U < u < 1.5p: 


7 7 f T { x 
NA kh akan Xx 
0.75 pabaya 
sheng: ah 
| \ 
salagah 
mae / + E + \ 
ose 27 cf a gi To t 
> LANE RON ATT eh RAE PA PAA pA df? a 
= ENNA ee eee eee oe ee 
e ee Che) Cee AG eee AA AA. 7 
DO KA KAM TAL AP A TTP AP PRE PE PE w/U + 1/2 
WA WA WA AA AA PEE ee Pr res 
NA WA MA UA AA 
oof NON EA WA! AI TG! Na! tiea! MI 
Yh ki TINA ET Tt ta Part oe ATE 
\ | | | | fi / 
075 UE 025 $ 025 05 075 
Re(z) 
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AJETET One-site Hubbard model 


Semiclassical partition function 


Using complex classical solutions Zm, let us calculate 


= 1 
Za = ` e7 Sm = e 50009, (« C EH 5) Aa ji 


m=— CO 


This expression is a good approximation for —1/2 < u/U 8 3/2. 
Nan 


[U 
-1/2 0 1 pg 
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AIETE One-site Hubbard model 


Sign problem after Lefschetz-thimble deformation 


This computation means that the sign problem exists after the 
Lefschetz-thimble deformation. 
Compute the reweighting factor: 


Is Tm dze 5") a" 03 (7 G oF 3) wa 


deje 6 Ep e TE 0 (0, 0270) 


eee 


At u = O, for example, 


Sop In ETO ao WA 
Fe z, ldz RS) T OP 8] 


Lefschetz-thimble method reduces Af = 2U to Af = U/8 in the 
one-site Hubbard model. 
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Maes GED 
Multi-flavor massless QED, 


2-dimensional U(1) gauge theories with Ny massless fermions: 


Ng 
Z= J DA e54 I DYDyexp | -X ii Pap, [Da - Ha] ba 
a=1 


Since the fermions are massless, the nonzero topological sectors do 
not appear: 

27 o, 27 1 

gt + z inde + *d@ + dr. 


o photon gauge 


A= 


o-dependence is computable using the anomaly equation, and does 
not have the sign problem. 


Yuya Tanizaki (RIKEN BNL) Fermionic sign problem Dec 12, 2016 @ IPMU 21 / 25 


ATE ntg Massless QED9 


Toron-field integral of multi-flavor massless QED» 
The toron-field integration becomes (7 = L/: temperature) 


27 
Z= i dhodh, exp |--Z Fto, m) A 
0 


where F is the fermion one-loop free energy (pw, = Lua/(27)), 


N 
1 ? T Lo 27 | 1 : 
F=N; (1-5) - EY fin (ear) 


621 n=l 


+ In (2 + gr AN) tn G 4 ee (nth —1+H4)4 aa) 
Hin G + er ria) ae) Ja 


In the limit 7 — 0, we can use the mean-field approximation with 
complex saddle points (to appear, with Motoi Tachibana) (cf.1504.02979, with 
Hiromichi Nishimura, Kouji Kashiwa). 
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Gradient flow for massless QED» 


All the relevant complex saddle points have real F: 


ZI NN 
HIN Hs 


bone 


Wz nl 
BK NI 
AMANN 
FF WA HABAN N NI 


Figure: Gradient flow in the Re(hi)- rIm(h 0) plane at Ng -3 


(to appear, with Motoi Tachibana) 


a a = = = Nao 


Phase structure of multi-flavor massless QED» 
At the zero-temperature with finite L, 7 = 0, the first-order 
transition occurs in this model: 
Ng] — Ng 


i bol 
27 
6 
4 
2 
0 31323 


/ 
7 "I 
2 3 3 
(a) Fermion number at Ny = 2 (b) Phase boundary at Ny = 3 
Since the mean-field approx. is good, one can show that 


2p .-5(h) 
Is. Aa Phe 
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Fermionic sign problem 


Summary 


ə The sign problem is reviewed from the viewpoint of exponential 
complexity. 


ə Using Cauchy's theorem, one can deform the oscillatory integral 
into the sum of steepest descent integrals on Lefschetz thimbles 


ə Monte Carlo method on Lefschetz thimble is developing. 
This has a potential to study not only the sign problem but also 
resurgence relation of QFTs numerically. 


e We consider two examples of the fermionic sign problem. 
In the one-site model, the exponential complexity is not solved, 
but Af is reduced. 
In massless QED», the exponential complexity is circumvented 
by using Lefschetz thimbles. 
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